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Abstract 

The radiation emitted by a charged particle moving along a helical orbit inside a dielec- 
tric cylinder immersed into a homogeneous medium is investigated. Expressions are derived 
for the electromagnetic potentials, electric and magnetic fields, and for the spectral-angular 
distribution of radiation in the exterior medium. It is shown that under the Cherenkov 
condition for dielectric permittivity of the cylinder and the velocity of the particle image 
on the cylinder surface, strong narrow peaks are present in the angular distribution for the 
number of radiated quanta. At these peaks the radiated energy exceeds the corresponding 
quantity for a homogeneous medium by some orders of magnitude. The results of numeri- 
cal calculations for the angular distribution of radiated quanta are presented and they are 
compared with the corresponding quantities for radiation in a homogeneous medium. The 
special case of relativistic charged particle motion along the direction of the cylinder axis 
with non-rclativistic transverse velocity (helical undulator) is considered in detail. Various 
regimes for the undulator parameter are discussed. It is shown that the presence of the 
cylinder can increase essentially the radiation intensity. 



PACS number(s): 41.60.Ap, 41.60. Bq 



1 Introduction 

The unique characteristics of synchrotron radiation, such as high intensity, high collimation, 
and the wide spectral range (see Refs. 03 El |IJ 03 E] and references therein), have resulted 
in extensive applications of synchrotron radiation in a wide variety of experiments and in many 
disciplines. These applications motivate the importance of investigations for various mechanisms 
of controlling the radiation parameters. From this point of view, it is of interest to consider the 
influence of a medium on the spectral and angular distributions of the synchrotron emission. 
This study is also important with respect to some astrophysical problems [HI El- The presence of 
medium can essentially change the characteristics of the high energy electromagnetic processes 
and gives rise to new types of phenomena. Well-known examples are Cherenkov, transition, 
and diffraction radiations. The operation of a number of devices assigned to production of 
electromagnetic radiation is based on the interaction of high-energy electrons with materials 
(see, for example, |Sj). 

*Email address: saharyan@server.physdep.r.am 
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The synchrotron radiation from a charged particle circulating in a homogeneous medium was 
considered by Tsytovich in Ref. (§] (see also Refs. |1UI lllj ). where it had been shown that the 
interference between the synchrotron and Cherenkov radiations leads to interesting effects. New 
interesting peculiarities arise in the case of inhomogeneous media. In particular, the interfaces of 
media can be used to control the radiation flow emitted by various systems. In a series of papers 
started in Refs. j!21 113j we have considered the most simple geometries of boundaries, namely, 
the boundaries with spherical and cylindrical symmetries. The synchrotron radiation from a 
charge rotating around a dielectric ball enclosed by a homogeneous medium is investigated in 
Refs. |13[ I14j. It has been shown that the interference between the synchrotron and Cherenkov 
radiations leads to interesting effects: if for the material of the ball and the particle velocity the 
Cherenkov condition is satisfied, there are strong narrow peaks in the radiation intensity. At 
these peaks the radiated energy exceeds the corresponding quantity in the case of a homogeneous 
medium by some orders of magnitude. A similar problem for the case of the cylindrical symmetry 
has been considered in Refs. ^21 EES! E3 El EES! In Ref. we have developed a recurrent 
scheme for constructing the Green function of the electromagnetic field for a medium consisting 
of an arbitrary number of coaxial cylindrical layers. The investigation of the radiation from a 
charged particle circulating around a dielectric cylinder immersed in a homogeneous medium, 
has shown that under the Cherenkov condition for the material of the cylinder and the velocity 
of the particle, there are narrow peaks in the angular distribution of the number of quanta 
emitted into the exterior space. For some values of the parameters the density of the number of 
quanta in these peaks exceeds the corresponding quantity for the radiation in vacuum by several 
orders. The radiation by a longitudinal charged oscillator moving with a constant drift velocity 
along the axis of a dielectric cylinder immersed in a homogeneous medium is investigated in 
Refs. |191 12(J| . It has been shown that the presence of the cylinder provides a possibility for the 
essential enhancement of the radiation intensity. 

In the present paper, on the basis of the Green function obtained in Ref. U2]> the radiation 
from a charged particle moving along a helical orbit inside a dielectric cylinder is studied. The 
corresponding problem for the charge moving in vacuum has been widely discussed in literature 
(see, e.g., Refs. El and references given therein). This type of electron motion is used 
in helical undulators for generating electromagnetic radiation in a narrow spectral interval at 
frequencies ranging from radio or millimeter waves to X-rays |22( I23j (see also El I24( I25j). 
The paper is organized as follows. In next section the vector potential and electromagnetic fields 
are determined for the region outside the cylinder by using Green function. A formula for the 
angular-frequency distribution for the radiation intensity is derived and various limiting cases are 
investigated in Sec. El In Sec. 0]we discuss the characteristic features of the radiation intensity 
and the results of the numerical evaluations are presented. Furthermore, we demonstrate the 
possibility for the appearance of strong narrow peaks in the angular distribution of the radiation 
intensity at a given harmonic. The case of helical undulator with non-relativistic transversal 
motion is considered in Sec. EJ Various regimes for the undulator parameter are investigated. 
Section El concludes the main results of the paper. 



2 Electromagnetic potentials and fields in the exterior region 

Consider a point charge q moving along the helical trajectory of radius po inside a dielectric 
cylinder with radius p\ and with permittivity eq. We will assume that this system is immersed in 
a homogeneous medium with dielectric permittivity e\ (for simplicity the magnetic permeability 
is taken to be unit). The particle velocities along the axis of the cylinder (drift velocity) and 
in the perpendicular plane we will denote by vn and v±, respectively. In a properly chosen 
cylindrical coordinate system (p, (ft, z) the corresponding motion is described by the coordinates 
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p = Po, 4> = u t, Z = V\\t, (1) 

where the z-axis coincides with the cylinder axis and ujq = vj_/po is the angular velocity of the 
charge. This type of motion can be produced by a uniform constant magnetic field directed 
along the axis of a cylinder, by a circularly polarized plane wave, or by a spatially periodic 
transverse magnetic field of constant absolute value and a direction that rotates as a function 
of the coordinate z. 

The solution to Maxwell's equations for the four-vector potential is expressed in terms of the 
Green function Gu(r,t,r',t') of the electromagnetic field being a second-rank tensor: 

Ai{v,t) = - JL_ J Guiv^v'^hiv'^dr'dt', (2) 

where the summation over / is understood and the components of the current density created 
by the charge are given by formula 

ji = ^vi5(p - po)5((f) - uj t)5(z - v\\t). (3) 

For static and cylindrically symmetric medium the Green function can be developed in a Fourier 
expansion: 

°° />oo 

Gn(r,t, r't') = / dk z duGii(m,k z ,u, p, p')exp[im((f)— (j)')+ik z (z — z') — iu(t — t')]. (4) 

The substitution of expressions (j2J) and (J3J into formula © gives the following result for the 
components of the vector potential: 



00 roo 



/U--X v (5) 

x [vxGty(m, k z ,u m (k z ),p,p ) + v\\Gi z (m, k z ,u m (k z ), p, p )] , 

where 

^m(kz) = mu + k z v\\. (6) 
Using the expression of the Green function given earlier in Ref. and introducing the notations 

X 3 = ~2 £ J- k zi J =0,1, (7) 

for the corresponding Fourier components G# = Gu(m, k z , u> m (k z ), p, po) in the region outside 
the cylinder, p > pi, we obtain 

G p4> = \ Y, pB ( % ] H m+p {\ lP ), (8a) 
p=±i 

ik z J m (\ po)H m (\ lPl ) s-^ J m+p (X pi) 

G pz = — > — — H m+p {Aip), (8b) 

2 P i a m W m f-£ Wm+p 



~ B$H m+p (\ lP ), (8c) 



2 p=±l 



= — > p — — H m+p (\ip), (8d) 

2pi a m Vv m Vy m+ p 

G« = ^^H m (X lP ), (8e) 

Pi 
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where J m (x) is the Bessel function and H m (x) = Hm\x) is the Hankel function of the first kind. 

(p) 

The coefficients Bm , p = ±1, in these formulae are determined by the expressions 



-'—'<rn 



1 



PiW, 



where 



m+p 



a, 



~ n "' l=±l 



rn+l 



SO 



E\ — £o 

Here and below we use the notation 



{=±1 



W, 



rn+l 



dH m (\ip x ) ,dJ m {\ QPl ) 

Jm\AQPl) E tlm{MPl) 



0) 



(10) 
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dpi dpi 

In the definition of Ai from Eq. (JJJ) one should take into account that in the presence of the 
imaginary part s"(uj) for dielectric permittivity (e\ = e'^ + ie") the radiation field in the exterior 
medium must damp exponentially for large p. This leads to the following relations: 



A I 



(u Wc)Vei - k^ c 2 /^, uj 2 m ei/c 2 > k 2 , 
W k z -u^i/c 2 , w^ei/c 2 < k 2 . 



(12) 



Note that for A? > the sign of Ai may be also determined from the principle of radiation 
(different signs of cot and Aip in the expressions for the fields) for large p. 

Substituting the expressions for the components of the Green function into formula we 
present the vector potential in the form of Fourier expansion 



°° /*oo 

Mr,t)= Y, e im ^-»ri dk z e lk ^- v ^A ml {m,k z ,p), 

m — ~, J— oo 



m=— oo 



where the Fourier components A m i = A m i(m,k z , p) are determined by formulae 

A r , 



qi 
2^ 



A,. 



2 P C$H m+p (\ iP ) 
p=±i 

p=±l 
q «[| 

ire piW r , 



-Jrn(^OPo)H m (Xlp), 



with the coefficients 



^m(AoPo)^m+p(Ao/5l)-ff m (AlPl) 



(13) 

(14a) 
(14b) 
(14c) 

(15) 



a cp±a m W m W m+ p 

The strengths of electromagnetic field are obtained by means of standard formulae of electrody- 
namics. As is seen from formula (|13|) . analogous expressions may also be written for the electric 
and magnetic fields. Having the vector potential, one can derive the corresponding Fourier 
coefficients for the magnetic field: 



nip 



iqk z 

2vr 



£ D$H m+p (A lP ), 



p=±i 



m</> 



H„ 



2ir 

p=±i 

2vr 



^ Y, pD$H m (\ iP ), 
p=±i 



(16a) 
(16b) 
(16c) 
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where the notation 



D (P) - C <P) - ^ Al Jm ( A °P°) n-±l 
ck z pi W m 



(17) 



is introduced. By making use of Maxwell's equation V x H = —iujei'E/c, one can derive the 
corresponding Fourier coefficients for the electric field: 



E, 



qc 



mp 



£ P H m+p (X lP ) \(^l + k A flW - a?dh») 



-En 



4vro; m (A; 2 )ei 



qic\\k z 



2TTuj rn (k z )ei 

v ' p=±i 



^ D$H m {\ x p). 



(18a) 
(18b) 
(18c) 



As follows from these formulae, E m • H m = 0, i.e., the Fourier components of the electric and 

magnetic fields are perpendicular. In the case of a particle motion in a homogeneous medium 

(p) (p) 

one has Eq = e\ and the terms in the definitions of the coefficients Cm and B m involving the 
function a m vanish. By taking into account that in this case W m = 2i/npi, one finds 



JT)(P) 



IT 

27 



v± 



Jm+p(^lPo) 



vn Ai 

ck 2 



-Jm(^lPo) 



(19) 



As functions of k z the Fourier coefficients for the fields determined by relations ()16[) . I|18[l 
have poles corresponding to the zeros of the function a m appearing in the denominators of 
Eqs. Pl).l|15|). It can be seen that this function has zeros only for A| < < Aq. As a necessary 
condition for this one has E\ < £q. Note that for the corresponding modes the Fourier coefficients 
are proportional to the MacDonald function K u (\\i\p), v = m, mil, and they are exponentially 
damped with the distance from the cylinder axis. These modes are precisely the eigenmodes 
of the dielectric cylinder and propagate inside the cylinder. Below, in the consideration of the 
intensity for the radiation to the exterior medium, we will neglect the contribution of the poles 
corresponding to these modes. 



3 Spectral-angular distribution of the radiation intensity 

Now we proceed to the consideration of the intensity for the radiation to the exterior medium. 
As it follows from the expressions of the fields, for Af < the corresponding Fourier components 
are exponentially damped for large values p, and the radiation is present only under the condition 
Xf > 0. The average energy flux per unit time through the cylindrical surface of radius p coaxial 
with the dielectric cylinder is given by the Poynting vector S: 

2-7T f T c 2tt 

I = — dt S • Tippdz, S = — E x H, T = — . (20) 
T J y_oo 4vr cj 

Substituting the corresponding Fourier expansions and using the fact that the replacement 
m — ► — m, k z — ► — k z leads to ui m (k z ) — ► —u> m (k z ), X± — > — Ai, H m (—k z ) = H^(k z ), we obtain 

oo „ 

/ = 2cnp Re £ ' / dk z [E m(j> H* mz - E mz H* m(j) ] , (21) 

m=0 ^ A ?>0 

where the prime over the sum means that the term with m = must be taken with the weight 
1/2. At large distances from the cylinder, by using the asymptotic expressions of the Hankel 
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functions, for the radiation intensity one finds 



a 2 r 2 ^ 
m=0 



dk. 



a2>o sx\u m (k z )\ 



-El 



£)(+!) _ £)("!) 



(22) 



where the coefficients Dm are defined by formulae (|17|). 

The case wo = for a fixed po corresponds to a charge moving with constant velocity vi\ on 
a straight line p = po parallel to the cylinder axis. For this case u> m (k z ) = k z v» and Eq. JJJ 



takes the form A 2 



1). Here and below we use the notations 



Pi 



-y/£j, (3j_ 



c 



(23) 



From the condition A| > it follows that the radiation is present only when the Cherenkov 
condition is satisfied for the parallel component of the particle velocity, vh, and dielectric per- 
mittivity Ei for the surrounding medium, where the radiation propagates: fan > 1. Introducing 
the angle $ of the wave vector with the cylinder axis, from the relation k z = u/vn we obtain 
cos i? = i-e., the radiation propagates under the Cherenkov angle of the external medium. 
As for the case under consideration v± = 0, the first term on the right of formula (|15|) vanishes 
and the dependence on po in the coefficients D$ is in the form of the Bessel function J m (\oPo)- 
It follows from here that in the limit po — > (particle moves along the axis of the cylinder) the 
term with m = contributes only and from Eq. Q22|) one obtains 



II 



Po=0 



4<7 v\\ 



J 



dw ■ 



1 



siUJ\W £ 



(24) 



where u> = \k 7 



\(0) 

W e = J (A( 0) pi)^i(Ai 0) pi) - ^Ji(Aj 0) P i)tfo(AfV), 

eiA 



(25) 



and 



A (o) 



(26) 



It can be checked that expression (|24|) coincides with the formula presented, for example in Ref. 
|26j . for the radiation of a charge moving along the axis of a cylindrical channel in a dielectric. 

Now we turn to the analysis of formula (1221) for the general case wo / 0. First of all we 
consider the contribution of the term with m = 0. For this term uj m (k z ) = k z Vi\ and, as in the 
previous case, from the condition Xf > it follows that the radiation at m = is present only 
when /%! > 1. This radiation propagates under the Cherenkov angle of the external medium. 

(p) 

By using the expressions for the coefficients Dm and introducing a new integration variable 
u = \k z \v\\, for the radiation intensity at m = one obtains 



4<? 2 V|| 

TT 2 pj 



duo 

n >l UE\ 



2 ^Jl{xfp,) J'^'po) 

Pl_|_ ,o h 



,(0). 



\w E 



(27) 



where Xj are defined by Eq. (J2 
substitute A 7 - = A^- . In the limit po 



and in the expression for W\ from Eq. we have to 

-> the first term in the square brackets vanishes and from 
Eq. (|27jl we obtain formula (|24jl for the radiation intensity from a charge uniformly moving along 
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the axis of a dielectric cylinder. In the case of the charge motion in a homogeneous medium one 
has Eq = e± and formula (|27|) for the radiation intensity at m = takes the form 



4°) = tf d ^ 

v \\ h n >i £ i 



/3iVi 2 (aS°Vo) + (^n - 1) J 2 (aS°Vo)] • (28) 



Next we consider the radiation intensity at m / harmonics. From the condition Af > 
one obtains the following quadratic inequality with respect to k z : 



,2, ,2 



- 13-, ) + 2k z ^ + ^ > 0. (29) 



Let the Cherenkov condition be not fulfilled initially for the drift velocity of the charge: /3]M < 1. 
In this case from inequality (|29|) we obtain: 



^V^TKir^Pn))- (30) 

It is convenient to introduce the new variable i? according to 

= B ^ (31) 

c 1 — pin cos if 

where from relation pUjl it follows that i? G (0,tt). Then from expression Q we have 

" wfe)= l^b ' - = 1,2... (32) 

Note that in accordance with Eqs. (|3T|) and (|32|) the quantities k z and u> m (k z ) are connected by 
the relation fc 2 = u) m (k z )y/£i cos i? /c. 

Now consider the case piy > 1, when the solution of inequality (|29|) has the following form: 



fc n^ , "3(SF«J u ("^rn)' 00 j- (33) 

Introducing again the variable $ according to expression (|3*Tj) . we see that $ G (0,i?o) U (^O) 71 ")) 
where i?q = arccos(/3^j| 1 ) is the corresponding Cherenkov angle for the drift velocity. The relation 
between the variables k z , u> m (k z ), ■& now is the same as in the case /3u\ < 1. At large distances 
from the charge trajectory the dependence of elementary waves on the space time coordinates 
has the form exp[uj m (k z )y/ei(psm$ + zcos# — ct / \/£i) / 'c], which describes the wave with the 
frequency 

uJ m = \u m {k g )\ = - T7, 34 

|1 — cosu| 

propagating at the angle i? to the z-axis. Formula (|34[) describes the normal Doppler effect in 
the cases p\ii < 1 and (3i\\ > 1, i? > $o, and anomalous Doppler effect in the case flu > 1, 
•& < i?o- By making use of the formulae given above, the expressions for A,- can be written as 

mw y/eo-ei cos 2 i? 
Ao = : 7. 5-, 35a 

c 1 — Pi|| COS V 

Ai = z 3 q • (35b 

c 1 — Pi\\ cos 17 
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Introducing instead of k z a new integration variable $ in accordance with Eq. (|31|) . from Eq. 
(I22fl for the radiation intensity atm/0 harmonics one finds 



E 

m=l 



dn 



dl„ 



dn 



(36) 



where dn = sin •dd'&dif) is the solid angle element, and 



dl m 
dfl 



,2, ,2„ 



2tt 3 c|1 - (3 n costf| 3 



D W_ D (-i) 



+ 



£,(1) + £,(-1) 



cos 2 1? 



(37) 



is the average power radiated by the charge at a single harmonic m into a unit solid angle. In 
the limit po ~^ for a fixed Cc>o the radiation intensity at the harmonic m ^ tends to zero as 
pQ m and the contribution from m = remains only. The latter corresponds to the Cherenkov 
radiation from a charge in uniform motion along the axis of the cylinder. 

In the case of a particle moving in a homogeneous medium with dielectric permittivity e±, 
in formula (|37|) one has £q = E\ and using formula (|19|) for the coefficients Dm > we obtain 



d/ (0) 

dn 



2fircy/ei\l - cosi?| 3 



AVm (AiPo) + 



cos # — 
sini? 



•^m(^lPo) 



(38) 



For ei = 1 this formula is derived in Ref. |2J (see also Refs. 0I3JI3]). In the case v\\ = and for 
an arbitrary e\, from (|38j) we obtain the formula for the synchrotron radiation in a homogeneous 
medium (see Refs. (HI El)- F° r a charge with a purely transversal motion (v» = 0) from Eqs. 

(fTH|). (fT7jl we see that L>^ } = (v/c)B$, and from Eq. (jSZI) the formula derived in Ref. dU is 
recovered. 

From the analysis presented above it follows that in presence of a dielectric cylinder the total 
intensity of the radiation can be presented as the sum 



I = h + l 



(39) 



where the first term on the right is given by formula (|27|) and describes the radiation with 
a continuous spectrum propagating at the Cherenkov angle of the external medium, if the 
condition [3i\\ > 1 is fulfilled (for (3^ < 1 this term is absent). The second term describes the 
radiation, which for a given angle $ has a discrete spectrum determined by formula (|34|) . With 
allowance for the dispersion of the dielectric permittivity, this term does not contribute to the 
radiation at the Cherenkov angle. This is connected with the fact that for a given m > the 
frequency defined by Eq. (|34|) . tends to infinity as $ approaches $0 and hence beginning from 
a certain frequency the Cherenkov condition ceases to be fulfilled. The angles for which the 
dispersion should be taken into account, are determined implicitly from the condition uo m > uo^ 
by using formula (|3"4"|) and frequency dependence of the permittivity e\ = £i(u m ), where oj<i is 
the characteristic frequency of the dispersion. 



4 Properties of the radiation 

In this section, on the base of formula from previous section we consider characteristic features 
of the radiation intensity. For a non-relativistic particle, /3i_l, f3\\\ <C 1, from the general formula 
one obtains 

dim _ 8g 2 c v /el / mp 1± \ 2{m+1) , 2 , 2{m _ x) 
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In this limit the energy in the higher harmonics is small compared to that in the fundamental, 
m = 1. 

For the case Eq < £ i the quantity Ao is purely imaginary in the angular region cos 2 i? > 
Eq/ei- By using the asymptotic formulae for the Bessel modified function for large values of 
the argument, it can be easily seen that in this region the radiation intensity exponentially 
decreases with increasing p\ in the limit when the wavelength for the radiation is much less than 
the cylinder radius: 

- 2 l A ^(|Ao|po), |Ao|pi » 1. (41) 



din 



dtt 



oc e 



This is caused by the fact that for Eq/ei < 1 the angle arccos(^/e /£i) corresponds to the angle of 
total internal reflection, and in the limit of the geometric optics the beams incident from inside 
on the cylinder surface cannot propagate at the angles cos 2 t9 > Eq/ei in the external medium. 

Now let us consider the behavior of the radiation intensity near the Cherenkov angle when 
|1 — cost?| <C 1. This behavior is radically different for the cases /3 || > 1 and /3 || < 1- In the 
first case the quantity Ao is real and using the asymptotic formulae for the cylinder functions 
for large values of the argument it can be seen that (assuming eq ^ e±) 



dl„ 



dQ 



OC |1 — p\|| COSt?| 4 , poll > 1. 



(42) 



In the second case, poll < 1, the quantity Ao is purely imaginary and the radiation intensity is 
exponentially suppressed: 



dim 

~dtt 



oc 1 1 — p\|i cos i?| 4 exp 



2muj (p 1 - po)^Jl - t 
1711 1 1 — f3u cost9| 



A)|| < i- 



(43) 



In the same limit, |1 — f3\\\ cost9| <C 1, for the intensity of radiation in a homogeneous medium 
from formula (|38|) one has 



dl {0) 

dn 



|1 - p n cost?| 



-2 



(44) 



Comparing with Eq. (|42|) . we see that in this limit the presence of the dielectric cylinder 
essentially increases the radiation intensity. Note that in the considered limit, in accordance 
with Eq. (|34|). the frequency for the radiation is large and the dispersion of the dielectric 
permittivity may be important. 

In Figure ^ we have plotted the dependence of the angular density for the number of the 
radiated quanta 

(45) 



dn 



hw m dQ. 



as a function on the angle < t? < ir for /3i± = 0.9, pi/po = 1.05, m = 10. The full and dashed 
curves correspond to the cases Eq/ei = 3 and £o/e\ = 1 (homogeneous medium) respectively. 
The left panel is plotted for p^y =0 (purely transversal motion, the curves are symmetric with 
respect to the rotation plane i9 = vr/2) and the right panel is for piy = 0.5. In both graphics the 
strong narrow peaks appear for the case when the dielectric cylinder is present. At the peaks of 
the left panel one has (Ty/Eihc/q 2 )dN m /dQ ~ 1.6- 10 and the width of the peaks is of the order 
At? ~ 6 • 10~ 5 . For the left peak on the right panel one has the value ~ 4.8 at the peak and 
At? « 0.01 and for the right peak one has the value 11 at the peak with the width At? 0.008. 



Figure |21 presents the angular dependence of the number of the radiated quanta for the values 
Pi/po = 1-5 (left panel) and pi/po = 3 (right panel). The values of the other parameters are 
p\_i_ = 0.9, m = 10, p^n = 0.5, Eq/ei = 3. For the left panel the value of the function at the peak 
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0.1 



0.08 



0.06 



0.04 
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Figure 1: The dependence of the angular density for the number of radiated quanta, 
(T^/eihe/q 2 )dN m /dQ, per period T of charge revolution as a function on the opening angle 
■& for = 0.9, pi/po = 1.05, m = 10. The left panel is plotted for (3\\\ = and the right panel 
is plotted for (3i\\ = 0.5. Full and dashed curves correspond to the cases £q/e\ = 3 and eo/ei = 1 
respectively. 




Figure 2: The dependence of the quantity (T 'y/elhc/ 'q 2 )dN m / 'd£l on the angle d for pi/ po = 1.5 
(left panel) and pi/po = 3 (right panel). The values of the other parameters are fii± = 0.9, 
m = 10, fan = 0.5, e /ei = 3. 
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corresponding to $ ~ 2.27 is equal ~ 2.5. For the left panel the peak corresponds to t? ~ 2.81 
with the value of the function ~ 11.5. 

To understand the presence of the strong narrow peaks in the angular distribution of the 
radiation intensity note that for large m from Debye's asymptotic expansion one has the following 
asymptotic formulae for the Bessel and Neumann functions 



with 



J m {my) ~ (signy) r 
Y m {my) ~ (signy) r 

C(y) = in 



exp [-m((y)] 
v / 2^(l-y 2 ) 1 /4' 

+1 2exp[mC(j/)] 
v / 2^(l - y 2 ) 1 / 4 



\y\ < i, 
\y\ < h 



(46) 



i + 



\y\ 



(47) 



For \y\ > 1 the corresponding formulae are obtained replacing the exponential functions by 
sin and cos (see, for instance, |2Z]). From formulae (|46|) it follows that for \y\ < 1 the ratio 
Jm{ m y) /Ym( m y) is exponentially small for large values m. Assuming that |Ai|pi < m and 
expanding over the small ratio J m {\\p{)/Y m (Xipi), for the coefficient a m defined by Eq. (jTU|) . 
one obtains 



Or. 



El — £q 



1 



i=±l 

AlJ m +z(AoPl)Jm(Al£>l) 



^m+;(Aopi) 
Y m+ i(Xipi) 



where 



AiJ m +z(Aopi)5^i(Aipi) — AoJ m (Ao/?i)^m+z(Aipi) 
Ai J m +i(Aopi)l" m (Aipi^ 



(48) 



F (0 

-*- TY1 



1 



(49) 



.Ao Jm(XoPl)Y m+ i(Xipi) 

and the second summand in the figure braces is exponentially small for large m. From here it 
follows that, at points where the real part of the function a m is equal to zero, the contribution 
of the imaginary part into the coefficients can be exponentially large. The corresponding 
condition for the real part to be zero has the form 



£4? 

i=±i 



2e 



E\ — Eq 



(50) 



Note that this equation is obtained from the equation determining the eigenmodes for the 
dielectric cylinder by the replacement H m — > Y m . 

For the further analysis it is convenient to rewrite equation (|5()|) in the form 



Ai 



Y' 

1 m 



(AiPi) 



A 



J; 



m(AoPl) A 



^m(AlPl) Jm{^0Pl)J V £®Y m {XiPi) 



A, 



ei l^(AiPi) 



A 



m 



A| 
A? 



^m(ApPl) 
■^m(AoPl) 

_ £i 



(51) 



First consider the case Aq < which is possible only for eq < £ l- Introducing instead of function 
Jm(Xopi) the Bessel modified function I m (\Xo\pi) and making use the corresponding uniform 
asymptotic expansions for large values of the order, we can solve equation ()51|) with respect to 
£o/ei for given Ao and Ai. As a result we obtain that e\ < Eq which is in contradiction with 
the condition Aq < 0. Hence, the above mentioned possibility for the appearance of peaks is not 
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realized for the case A 2 , < which corresponds to the angular region cos 2 ?? > £q/£\. For this 
reason we will concentrate on the case cos 2 $ < Sq/si. 
Let us consider separate cases. 

a) Let firstly (AoPo/ m ) 2 > 1 and, hence, (Xopi/m) 2 > 1. By taking into account the condition 
|Ai|pi < m, assumed earlier, we obtain that Af < Aq and, hence, £o > £\. For the solutions to 
equation (|5U[) the contribution from the first summand in the figure braces of Eq. (|48[) vanishes. 
By taking into account asymptotics (|46[). for the angular density of the number of the radiated 
quanta we obtain the following estimate 



dN m 



oc exp [2m£(Aipi/m)] , (52) 



with the function £(y) from Eq. (|47[). As this function is monotonic decreasing the exponent 
for a given m is larger for smaller values of |Ai|. 

b) Let now (Xopo/m) 2 < 1 and (Xopi/m) 2 > 1. The substitution of the asymptotic formulae 
(jH^ into Eq. (j37|) shows that, the peaks arise only under the condition (Xipi) 2 < (AoPo) 2 ; which 
is possible only for eo > £\. For the angular density of the radiated quanta at the peaks one 
obtains ^ 

oc exp {2m {((X 1 p 1 /m) - ((X p /m)}} . (53) 

c) In the case (Xopi/m) 2 < 1 we have also (Xopo/m) 2 < 1. Again, using the asymptotic 
formulae (|4fi(l. it can be seen that in this case Eq. (|51(1 has no solutions and the peaks are 
absent. 

Therefore, as necessary conditions for the presence of the strong narrow peaks in the angular 
distribution for the radiation intensity one has |Ai|/?i < m < |Ao|pi, or in terms of the angle $ 



l ^JT l sm$< II -All costf| < ^Veo - £i cos 2 . (54) 
c 11 c 

From these conditions it follows that we should have 

£ >£i, vVeo/c> 1, (55) 

where v = ^v 2 + uj^p 2 is the velocity of the charge image on the cylinder surface. In particular, 

the second condition in Eq. (|55|) tells that the Cherenkov condition should be satisfied for the 
velocity of the charge image on the cylinder surface and dielectric permittivity of the cylinder. 
When the Cherenkov condition for the velocity of the charge image and dielectric permittivity 
of the surrounding medium is not satisfied, v^/Fi/c < 1, the possible strong peaks are located 
in the angular regions defined by the inequality 



— yj£\ COS 17 

C V 



PlV± V 2 

< — -V^o-l- ( 56 ) 
pov V c z 



If Cherenkov condition Vy/el/c > 1 is satisfied, in addition to this inequality the condition 



-v/el cos $ 

c v 



> ea±je Sl - 1 (w) 

pov V c z 



has to be satisfied. By the presented arguments we can also estimate the width for the peaks. 
Taking into account (flS]) and expanding a m near the angle corresponding to the peak, ■& = $ p , 
it can be seen that the angular dependence of the radiation intensity near the peak has the form 

dim, 1 / dim, i 

(58J 



dn (# - $ p ) 2 /b 2 p + i\ dn y „ , 
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where 

b p ~ J m +i(^iPi)/y m +i(^iPi)\i)={) p ~ exp[-2mC(Aipi/m)]. (59) 
It follows from here that the width of the peak has an order 

A#~ exp[-2mC(Aipi/m)]. (60) 

As we see from estimates (j52j) . (|53|) . at the peaks the angular density for the number of 
emitted quanta exponentially increases with increasing m. However, one has to take into ac- 
count that in the realistic situations the growth of the radiation intensity is limited by several 
factors. In particular, the factor which limits the increase, is the imaginary part for the dielectric 
permittivity e" , j = 0, 1. This leads to an additional summand in the denominator of formula 
(|58|) , proportional to the ratio e" / e'j , where e'j is the real part of the dielectric permittivity. As 
a result for bp < £j/£j the intensity and the width of the peak is determined by this summand 
and the saturation of the radiation intensity at the peak takes place. 



5 Helical undulator 

Let us consider an important case relativistic charge motion in direction of the cylinder axis 
when the velocity of orthogonal motion is non-relativistic, j3i± <C 1. In helical undulators this 
type of motion is realized by a spatially periodic transverse magnetic field of constant absolute 
value and a direction that rotates as a function of the longitudinal coordinate z. In the discussion 
below we will assume that < 1 — <C 1. Under this condition, as it follows from Eq. ()37|) . the 

radiation intensity is sharply peaked for opening angles near ~ l/7i|| , where 7^ = l/^/l — (3^ 

is the longitudinal Lorentz factor. It is convenient to introduce new variables 

X = 7i||tf> « = /3l±7i||> ( 61 ) 
in terms of which for small angles 1? one has 



1 + x 2 .„ X 



1 - /^H cos •& » n _ 2 A , dft = \dxd<\>. (62) 



The parameter k is the so-called undulator parameter which characterizes the relation between 
velocity deviation angle and the opening of the radiation cone. The cases k > 1 and k < 1 
correspond to wigglers and undulators, respectively (strong and weak undulators in another 
terminology). For the quantities X1P0 and X0P0 we obtain 

2rriK,x 2m«7 1 || , — 

1 + X 1 + X 

On the base of these relations three regimes may be distinguished. 

a) The first one corresponds to the case k > 1 when \\po > 1 and A0P0 ^ 1- I n the defini- 
tion of the coefficients D$ the Bessel functions J^(XoPj) can be replaced by their asymptotic 

expansions for large values of the argument (we assume that po an d pi are of the same order). 

(p) 

The main contribution into the coefficient Dm comes from the second summand on the right 
of formula (|15|) . This term dominates the contribution coming from other terms in Eqs. (|15jl. 
(|17|) by the factor 7111. It can be seen that the radiation intensity as a function on $ has sharp 
peaks at the points $ = $1 corresponding to the zeros of the Bessel function J m (\oPi), i-e. at 
the angles defined from Xopi = j m i, where j m ,z is the Zth positive zero of the function Jm{z). 
Note that as Ao/>i 3> 1, for the zeros j m j we can use McMahon's expansion for large zeros (see, 
for instance, [57]). For the radiation intensity at the corresponding peaks one has an estimate 
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dli / d$d(j) oc 7^1, . At the angles away from the peaks the radiation intensity is relatively sup- 
pressed by the factor 7x11: dli/dfidcj) oc 7^,,. This is illustrated in Figure El where the number of 
radiated quanta at harmonics m = 1 and m = 2 is presented as a function on the variable \ f° r 
the value of the undulator parameter k = 1. We verified numerically that the sharp peaks in 
these graphs are located at the angles corresponding to the zeros of the function J m (Ao/5i). In 
this regime the particle is a relativistic one in the reference frame moving with the velocity v\\ 
and the spectrum contains many harmonics. 




350 ■ 




Figure 3: The dependence of the quantity (T y /eihc/q 2 )dN m /d'&d(j) on the reduced angle % f° r 
m = 1 (left panel) and m = 2 (right panel). The values of the other parameters are 7111 = 10, 
Pi/po = 1-5, k = 1, eq/ex = 3. 



b) In the second regime one has ivy\n ~ 1 which corresponds Aipo ^ 1; A0P0 ~ 1- I n this case 
we replace the Hankel functions in the definitions of the coefficients by their asymptotic 
expressions for small values of the argument. This leads to the estimate dli/d'ddcj) oc 7i||/ In 2 7111 
and dlm/d'dd^ oc l/7 2 jj" _3 for the harmonics m > 2. As we see in this case the main contribution 
comes from the fundamental harmonic m = 1. This is a simple consequence of the fact that 
under the condition k <C 1, the particle has a non-relativistic velocity in the reference frame 
which moves with the velocity vu. 

c) The third regime corresponds to k <C 7^ji when both quantities Xipo and Ao/>o are small. 
Expanding the cylindrical functions, we can see that, as in the previous case, the main part of 
the radiation occurs at the m = 1 harmonic and to the leading order one receives 



dh 



2 c(«7i 



2 3/2 



£q — E\ 



X 



Eo + £lJ (1 + X 



2\5 



dxd(f>. 



(64) 



For the higher harmonics one has an estimate dlm/d'dd^ oc («7i||) l^vt ' m — ^" I n t e g ra t 
ing over x and we find the total power of radiation in this regime: 



2q 2 c 

2 3/2 

Po e i 



(«7l||) 4 



£0 + £i 



(65) 



As in the case a), in regimes b) and c) the main contribution into the radiation intensity comes 
from the second term on the right of formula (j!5j) . 

For a particle moving in a homogeneous medium, £0 = E\, the leading term vanishes and it 
is necessary to take the next terms in the corresponding asymptotic expansion. For small values 
of the undulator parameter, k -C 1, the radiation intensity is dominated by the basic harmonic 
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and the corresponding angular distribution is given by the formula (see, for example, Ref. [S] 
for the case £\ = 1) 

(0) _ 2q 2 CK 4 1 + x 4 



dl\ 



Xdxdcp- 



(66) 



Integrating over the angles we obtain the total power of the radiation in a homogeneous medium: 



AO) „ 2 1 c 

1 ~ o 2 3/2' 

Now comparing formulae l|65|) and (|67|) . we see that 



(67) 



- 3 7i|| 



£0 + £l 



(68) 



and the presence of the cylinder essentially increases the radiated power to compared with the 
radiation in the homogeneous medium. 

The results of the numerical evaluations of the radiation intensity for small values of the 
undulator parameter are presented in Figure where the function (T yfeihc/ q 2 )dN m / dfidcf) is 
plotted versus the undulator parameter k and the reduced angle x f° r the values of parameters 
7x|i = 10, m = 1, pi/po = 1.05. The left panel corresponds to the radiation from a charge 
moving in a homogeneous medium (eq = £%) and the right panel corresponds to the presence of 
the cylinder with dielectric permittivity Eq = 3ei. 





Figure 4: The dependence of the quantity (T y /Eihc/q 2 )dN rn /d'dd(j) on the undulator parameter 
k and the reduced angle x f° r £ o = £i (left panel) and £o = 3ei (right panel). The values of the 
other parameters are = 10, m = 1, pi/po = 1.05. 



6 Conclusion 

Synchrotron radiation has become one of the most valuable and useful scientific tools with ever 
increasing applications for basic and applied research. In the present paper we have considered 
the influence of a dielectric cylinder on the synchrotron radiation from a charged particle moving 
along a helical orbit inside a dielectric cylinder. By using the Green function, we have derived 
formulae for the electromagnetic potentials and fields, Eqs. I|14j). (|16|) . (|18|) . and for the angular- 
frequency distribution of the radiation intensity emitted into the exterior medium. The latter 
is a sum of two terms. The first one, given by formula (|27j). is present when the Cherenkov 



15 



condition for dielectric permittivity of the exterior medium and drift velocity of the charge is 
satisfied, > 1, and describes the radiation with a continuous spectrum propagating under 
the Cherenkov angle of the external medium. The second term in the total radiation intensity 
is given by formula ()36|) . where the separate terms (|37|) describe the radiation, which for a given 
propagation direction has a discrete spectrum determined by formula (|34|) . We have investigated 
various limiting cases of the general formula for the radiation intensity and have compared this 
formula to previously known special cases (radiation in vacuum, purely transversal synchrotron 
radiation, radiation from a charge moving with constant velocity along the axis parallel to the 
cylinder axis). In Section 0] we have investigated characteristic features of the radiation. Our 
numerical calculations have shown that under certain conditions strong narrow peaks appear 
in the angular distribution of the radiation intensity. By using asymptotic formulae for the 
cylindrical functions, we analytically argued the possibility for the presence of such peaks and 
the corresponding necessary conditions are specified. In particular, the peaks are present only 
when dielectric permittivity of the cylinder is greater than the permittivity for the surrounding 
medium and the Cherenkov condition is satisfied for the velocity of charge image on the cylinder 
surface and the dielectric permittivity of the cylinder. We have also estimated the width of 
these peaks. The peaks correspond to the opening angles for which the real part of the function 
a m from IjlOjl vanishes. The corresponding condition is in form (|5fl|) and is obtained from the 
equation determining the eigenmodes for the dielectric cylinder replacing the Hankel functions 
by the Neumann ones. In section 00 an application of the general formula is considered to 
the practically important case of the non-relativistic transverse motion when the longitudinal 
motion is relativistic. This type of motion can be produced by a variety of electromagnetic field 
structures and is used in helical undulators for generating electromagnetic radiation in a narrow 
spectral interval. Different regimes for the undulator parameter are considered corresponding to 
strong and weak undulators. In the former case the radiation intensity contains many harmonics 
and has sharp peaks at the opening angles determined by the zeros of the Bessel functions 
Jm(^oPi)- At the angles away from these peaks the radiation intensity is relatively suppressed 
by the factor 7^1 . In the weak undulator regime the radiation is dominated by the fundamental 
harmonic and the leading contribution into the intensity comes from second term on the right 
of Eq. (|15|) . In the case of a particle moving in a homogeneous medium this term vanishes and 
it is necessary to take the subleading terms. The corresponding numerical results are plotted in 
Figures |3] and These results show that the presence of the cylinder provides a possibility for 
an essential enhancement of the radiated power as compared to the radiation in a homogeneous 
medium. 
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